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ABSTRACT 


Detailed experimental investigation of second generation controlled-diffusion 
compressor stator blades at design inlet-flow angle was performed in a low-speed cascade 
wind tunnel using various experimental methods. Surface pressure measurements were 
obtained using three instrumented blades, from which coefficients of pressure were 
calculated. Laser-Doppler velocimetry was used to characterize the flow in the inlet, in 
the passage between two blades, in the boundary layer of the blades, and in the wake. A 
five-hole pressure probe was used to determine the loss coefficient and the axial-velocity- 
density ratio of the flow through the cascade. Although the blades produced significant 
lift, separated flow was discovered on the suction side of the blades at approximately fifty 
percent axial chord, which showed that the design was not totally successful. All the 
experimental measurements were performed at an inlet flow Mach number of 0.22 and a 
Reynolds number, based on chord length, of 640,000. 

Experimental blade-surface pressure coefficients were compared with values 
predicted using a computational fluid dynamics code. These initial predictions did not 
match well with the experimental results. 
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I. INTRODUCTION 


A. BACKGROUND 

Increased performance requirements for aircraft gas turbine engines necessitate 
continuous progress in gas turbine engine research and development. Engine compressor 
stall and off-design behavior continue to limit the performance of military aircraft gas 
turbine engines. Compressor stall can lead to degradation in engine performance or even 
total loss of engine power, which could result in mission failure or aircraft loss. 
Additionally, the increased thrust-to-weight ratio engines needed for present and future 
military and civil aircraft require improved compressor design for increased performance 
at the same or reduced weight. 

Compressor blading has progressed over the years from the use of NACA-65 
series. Double Circular Arc (DCA) and Multiple Circular Arc (MCA) bladeshapes to the 
design of Controlled-Diffusion (CD) blading. CD blading was developed to control 
diffusion on the suction side of the airfoil, thus avoiding boundary-layer separation at the 
engine design point. Transonic design allowed shock-free operation in the transonic range 
which increased compressor and engine efficiency. The resulting increase in efficiency 
could be used to improve engine performance, or the same engine performance could be 
achieved using fewer compressor blades, allowing a reduction in engine weight. 

The CD compressor blades investigated in the present study were designed by 
Thomas F. Gelder of NASA Lewis Research Center (LeRC) [Ref. 1]. The compressor 
stator profiles were Stator 67B blades, which together with Rotor 67 comprised 
Compressor Stage 67B. Compressor Stage 67B was previously studied experimentally at 
the NASA LeRC compressor test facility [Ref. 1], The Stator 67B blades were second 
generation CD blades which were designed as an improvement to the Stator 67A CD 
blades. The Stator 67A blades, designed by Nelson Sanger [Ref. 2], coupled with Rotor 
67 to form Stage 67A. Ten midspan Stator 67B blade profiles were machined from 
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aluminum and installed in the Naval Postgraduate School (NPS) Turbopropulsion 
Laboratory (TPL) Low-Speed Cascade Wind Tunnel (LSCWT) for the present study. 
Pneumatic probe measurements were made upstream and downstream of the blading. 
Laser-Doppler velocimetry (LDV) measurements were performed upstream of the blade 
row, in the passage between the blades, in the boundary layers, and downstream in the 
wake, at the design inlet-flow angle. A quasi-three dimensional computational fluid 
dynamics (CFD) code. Rotor Viscous Code Quasi-3-D (RVCQ3D) Version 300, was also 
used for comparison with the experimental results. 

The Stator 67A compressor blading was the focus of a decade of research at the 
NPS TPL. Koyuncu [Ref. 3] studied the performance at various incidence angles using a 
five-hole probe. Dreon [Ref. 4] improved the accuracy in loss and pressure measurement 
near design conditions. Elazar [Ref. 5] performed LDV measurements in the boundary 
layer, in the blade passage and in the wake at three incidence angles. Murray [Ref. 6] 
upgraded the LDV system and measured the flow in the wake near stall. Classick [Ref. 7] 
upgraded the pressure probe data acquisition software and performed surveys near stall. 
Armstrong [Ref. 8] further improved the data acquisition system and compared flow at 
high incidence angles over a blade with a modified leading edge and an unmodified blade. 
Hobson and Shreeve [Ref. 9] used LDV to study the flow at very high incidence. Ganiam 
[Ref. 10] performed LDV measurements at stall. Williams [Ref. 11] verified the LDV 
measurements at stall and obtained a prediction using computational fluid dynamics 
(CFD). 

B. PURPOSE 

The objective of the present study was to install the Stator 67B midspan blade 
sections in the NPS cascade and to examine the flow through the blading at the design 
inlet-flow angle using LDV and pressure probe measurements. The intent was first to 
obtain a thorough understanding of the flow at the design inlet-flow angle, and to validate 
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the design itself. Future studies of Stator 67B cascade blading at higher inlet-flow angles 
will determine the off-design and stalling behavior. 
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II. TEST FACILITY AND INSTRUMENTATION 


A. LOW-SPEED CASCADE WIND TUNNEL 

The NPS Turbopropulsion Laboratory Low-Speed Cascade Wind Tunnel Building 
is shown in Figure 1. Elazar [Ref. 5] thoroughly documented the uniformity of tunnel 
flow conditions and the periodicity of the flow in the cascade test section with 20 Stator 
67A blades in the cascade at 40 (design), 43 and 46 degrees inlet-flow angle. 


B. TEST SECTION 

A schematic of the cascade is shown in Figure 2. Ten CD blades with elliptical 
leading- and trailing-edges were scaled from the midspan section of Stator 67B [Ref. 1] 
and machined from aluminum using numerically-controlled machining methods. Figure 3 
shows the profile of the blade. Table 1 contains the machine coordinates used to 
manufacture the blades. Three blades were instrumented with pressure taps which were 
machined into fine metal tubing laid into the blade surface. A fully-instrumented blade 
containing 42 pressure taps was installed at blade location number 6 [Figure 2], Two 
blades numbered 2 and 8 were partially instrumented with eight pressure taps each. Figure 
4 shows the pressure tap locations on the instrumented blades. Additionally, blades 3 and 
4 were treated with a black anodized coating to minimize laser light scatter for the LDV 
measurements. The LDV test passage between blades 3 and 4 is shown in Figure 5 with 
the locations of the 13 measurement stations shown as fractions of the axial chord (C ac ). 
Table 2 contains a summary listing of the geometrical parameters of the cascade test 
section. 

The ten blades were installed in the test section using brass shims for alignment 
perpendicular to the cascade endwalls. A digital inclinometer with an accuracy of 0.1 
degrees was used to set the stagger angle of the blades. The tunnel inlet sidewalls were 
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adjusted to the design inlet-flow angle, Pi = 36.3 degrees. The inlet guide vanes were 
adjusted until the mean inlet flow to the test section at Station 1 was equivalent to the 
design inlet-flow angle. The final tunnel adjustment was that of the tail board setting 
angles which were adjusted to give uniform wall static pressures downstream of the 
blading [Figure 2], 


C. INSTRUMENTATION 

1. Pneumatic Data Acquisition System 

Two different pneumatic data acquisition systems were used with a traversing five- 
hole probe to determine the loss coefficient (co), the axial-velocity-density ratio (AVDR), 
and blade surface pressure coefficients (Cp). The first system was a computer-controlled 
automated data acquisition system and software documented by Classick [Ref. 7] and 
modified by Armstrong [Ref. 8], The Hewlett-Packard data acquisition system hardware 
and software “ACQUIRE” and “LOSS” are fully described by Armstrong [Ref. 8], 

The second data acquisition system consisted of a Scanivalve Digital Interface Unit 
(SDIU Mk5), a Scanivalve Controller (CTLR2/S2-S6), a Hewlett-Packard digital 
voltmeter (HP 343 7A) and a 48 channel Scanivalve rotary pressure scanner. Loss and 
AVDR calculations were performed using a personal computer with the spreadsheet 
software Microsoft “Excel”. The five-hole probe and the five-hole probe traverse used for 
both systems are shown in the photographs in Figures 6 and 7. Figure 8 shows a 
schematic of the second data acquisition system. Plenum pressure and temperature 
sensors were the same as those used in the first system. A Prandtl probe located upstream 
of the test section was used as the total pressure reference for the second system. 

Prior to implementation of either of the two acquisition systems, a series of surface 
pressure measurements was performed using banks of water manometers. The partially- 
instrumented blades, mounted at blade position numbers 2 and 8 in the test section, were 
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connected to manometers for an initial observation of the pressure distribution over the 
blades. 


2. Laser Doppler Velocimeter 

A four beam, two color TSI Model 9100-7 LDV system was used for all surveys. 
Elazar [Ref. 5] thoroughly described the LDV system, including the laser model, optics, 
atomizer and seeding, and data acquisition. Murray [Ref. 6] described the fully automated 
traverse mechanism. The LDV and traverse systems were controlled by a personal 
computer. Experimental data were processed using TSI FIND software on the personal 
computer. The LDV and traverse systems are shown in Figure 9. 
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Figure 5. Survey Station Numbering and Position 


Blade Type 

Stator 67B Controlled-Diffusion 

Number of Blades 

10 

Blade Spacing 

152.4 mm 

Chord 

127.14 mm 

Solidity 

0.834 

Thickness/Chord 

0.05 

Setting Angle 

16.3° ± 0.1° 

Span 

254.0 mm 


Table 2. Test Section Data. 
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Figure 7. Five-Hole Pressure Probe Traverse Mechanism. 
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III. EXPERIMENTAL PROCEDURE 


A. SURFACE PRESSURE MEASUREMENTS 

The tunnel was started and allowed to reach an equilibrium operating temperature 
prior to all surveys. Once the tunnel reached operating temperature, the plenum chamber 
pressure was set and maintained at 12 inches of water. Twelve inches of water for the 
plenum pressure gave an average inlet Mach number of 0.22 and an average Reynolds 
number based on chord, Re c , of approximately 640,000. 

1. Water Manometer System 

Prior to the automated data acquisition system experiments, surface pressure 
measurements were obtained using banks of water manometers connected to the two 
partially-instrumented blades. The water manometers were manually recorded to an 
accuracy of 0.1 inches of water. Plenum total pressure and temperature, wall static 
pressure, Prandtl probe total and static pressures, and local atmospheric pressure were 
recorded in addition to the blade pressure taps. Calculations of the coefficients of 
pressure, Cp, based on the Prandtl probe total and static pressures, were performed using 
the software Microsoft “Excel”. 

2. HP Automated Data Acquisition System 

The HP data acquisition system was initialized, and the Scanivalve was calibrated 
using a water manometer prior to each experimental run. A U-tube water manometer and 
a regulated compressed-air supply were used as the calibration reference for both of the 
Scanivalve rotary pressure scanners. The five-hole probe was positioned away from the 
three instrumented blades to prevent wake interference. The program “ACQUIRE” 
sampled and recorded pressures using a Scanivalve rotary pressure scanner, and then 
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calculated Cp for each blade pressure tap based on the plenum total pressure, as described 
by Armstrong [Ref. 8], 


B. FIVE-HOLE PRESSURE PROBE MEASUREMENTS 

The five-hole probe traverse mechanism was mounted upstream of the test section 
at the tunnel lower survey location for the loss and AVDR measurements. The five-hole 
probe was then mounted on the traverse mechanism and centered at the blade midspan 
position in the tunnel. A U-tube water manometer was used with a regulated compressed- 
air supply to calibrate the Scanivalve rotary pressure scanner prior to making the 
measurements. The five-hole probe yaw transducer was also calibrated using a digital 
inclinometer. The five-hole probe was positioned in the pitchwise direction and aligned 
with the tunnel flow by balancing the probe in yaw using a U-tube water manometer prior 
to taking pressure data. The five-hole probe calibration coefficients developed by 
Armstrong [Ref. 8] were used by the data analysis. 

1. HP Automated Acquisition System 

The HP automated acquisition system was used to determine loss and AVDR for 
the blade section as outlined by Armstrong [Ref. 8], The five-hole probe was positioned 
in both the lower and the upper traverse positions. Data were taken over a full passage 
width of 154 mm in increments of 2.54 mm, 5.08 mm and 12.77 mm, respectively, for the 
surveys using Armstrong’s “ACQUIRE” program [Ref. 8], The “LOSS” program was 
then used to calculate loss and AVDR for each combination of upstream and downstream 
data, as described by Armstrong [Ref. 8], 

2. Manual Loss Coefficient and AVDR Measurements 

The SDIU, Scanivalve controller, Scanivalve rotary pressure scanner and HP 
digital voltmeter were used with the five-hole probe to determine the loss coefficient and 
AVDR. The SDIU and controller were used to step the Scanivalve to each port for 
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measurements, and the pressure data were read from the digital voltmeter and recorded 
into a spreadsheet on a personal computer for data reduction. Plenum total pressure and 
temperature, atmospheric pressure, wall static pressure, and Prandtl probe total and static 
pressures were recorded. Probe yaw angle was also recorded. The spacing interval for 
the survey was 2.54 mm over a total survey width of 152.4 mm with a total of 61 data 
positions being recorded. 

Appendix A contains the formulas which were used for the loss calculations. The 
calculation method was as follows. First, the dimensionless velocity coefficients, 3, were 
determined from the five-hole probe pressure data. Next, the dimensionless velocities, X, 
were determined using the polynomial coefficients for the five-hole probe. The 
dimensionless reference velocity, X re f, was determined using the pressures from the Prandtl 
probe. The reference flow (massflux) functions, K, were calculated for both the upstream 
and the downstream data. AVDR was determined using a Romberg numerical integration 
routine from the values of K. The loss coefficient, to, was determined from values of C P K, 
using the same type of numerical integration scheme. 


C. LDV MEASUREMENTS 

1. LDV Probe Volume Alignment and Reference 

The need for a precise reference for the LDV probe volume necessitated the design 
of an alignment tool. A photograph of the tool, which was manufactured out of 
aluminum, is shown in Figure 10, and the tool machine drawing with dimensions is 
presented in Appendix B. The tool was positioned on blades 3 and 4 in the spanwise 
direction by placing the two end pieces against the back of the test section wall. The two 
endpieces secured the tool tightly against the trailing edge of both blades using the 
machine screw on the right-hand side, which ensured consistent positioning. Three laser 
alignment holes were machined in the tool’s center block, each with a 0.3302 mm 
diameter. The laser probe volume was focused through one of these small holes, and the 
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automated traverse table (TSI Model 9500) was initialized to the coordinates of the hole. 
This alignment process was performed for each survey. 

The zero reference point was determined using the blade geometry at the design 
setting angle. The reference point was defined by the normal intersection of two lines 
which were tangent to the minimum x and y points of the leading-edge ellipse of blade 
number 3. At this intersection, x and y were defined to be zero. The spanwise mid-point 
of the blade, at 127 mm, was used as the third zero reference point, and all measurements 
were performed at midspan. 

2. Inlet Guide Vanes Adjustment 

Initial LDV inlet surveys were performed over a 254 mm distance, or nearly two 
blade passage widths, at Station 1 to allow proper adjustment of the inlet guide vanes. 
The mean inlet flow angle was computed from survey data and the inlet guide vanes were 
adjusted until this angle converged to within 0.1 degree of the design angle of 36.3 
degrees. 

3. LDV Surveys 

A total of 28 surveys were completed, including at least one survey of each of the 
thirteen stations and four boundary layer surveys. All surveys were taken by collecting 
1000 data points at each location across the station, with the exception of the Station 1 
surveys, where 3000 data points were taken. Figure 11 shows the location of the station 
surveys which will be discussed in the body of this report. A listing of all experimental 
runs performed is included in Appendix C. 

The laser optics were configured identically throughout this study, with the 514.5 
nm green beam measuring the axial (or vertical) velocity component, U, and with the 488 
nm blue beam measuring the tangential (or horizontal) velocity component, V. Frequency 
shifting of 5 MHz was performed, as outlined by Ganaim [Ref. 10], to detect reverse flow 
velocities. Alignment was performed prior to each measurement using the laser alignment 
tool, as described above. 
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Ambient pressure, and plenum total pressure and plenum total temperature were 
recorded for each survey. The velocity at the inlet to the test section was used as the 
tunnel reference velocity, V re f. The reference velocity was used to nondimensionalize the 
velocity measurements. A series of calibration surveys at six plenum pressures and 
temperatures were performed and a least-squares curve fit was applied to the data to 
determine a calibration curve for Vref. A Newton-method iteration algorithm was used to 
determine V re r for each survey using plenum total pressure, plenum total temperature and 
atmospheric pressure. A FORTRAN code “CALIBI” was used to calculate V ref , and the 
FORTRAN code listing and input and output data files are contained in Appendix D. 

TSI Flow Information Display (FIND) Version 4.0 software was used to acquire 
and analyze all LDV data. Velocity, turbulence intensity, Reynolds stress and Reynolds 
stress correlation coefficient information were processed using FIND, and then further 
non-dimensionalized using the inlet flow reference velocity, V re f, to allow comparison with 
data taken at different LDV configurations (i.e. LDV axis not perpendicular to the optical 
access window). 

a. Inlet Surveys 

The inlet flow region was surveyed across Stations 1, 2 and 3, in the far- 
and near-upstream regions. The laser system was positioned horizontally for the Station 1 
and 2 measurements. Station 3 measurements were conducted first with no pitch, for a 
region comprising 93 percent of passage width, and then the laser was pitched upward 
four degrees, allowing the full passage to be studied. The potential problems due to 
pitching and yawing of the LDV system were discussed by Hobson and Shreeve [Ref. 9], 
Their analysis showed a maximum spatial error from probe volume orientation to be 0.3 
mm, which is the measurement volume minimum diameter. 

b. Passage Sundays 

Passage surveys were conducted at Stations 4 through 10. The laser was 
horizontal for all passage surveys. One survey taken at Station 7 was performed with a 
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laser yaw angle of 4 degrees, but no pitch. Initial and final measurement boundaries were 
determined for the passage surveys by experimenting with the laser position until an 
adequate data rate was achieved. Interference degraded the laser beams when positioned 
too close to the blade surface, thus diminishing the data rate of the LDV system to 
unusably low levels. The boundaries of the surveys were set to exclude regions of 
interference. 


c. Boundary Layer Surveys 

Boundary layer surveys were performed at one station on the pressure side 
and three stations on the suction side of the blade. The boundary layer surveys on the 
suction surface of blade 3 were performed at Stations 5, 7 and 9. The boundary layer 
survey on the pressure surface of blade 4 was performed at Station 8. Prior to the 
experiment, the laser alignment tool was used to align the laser at the measurement pitch 
and yaw angles. The laser pitch and yaw angles used for the surveys were determined 
experimentally by finding the angle which allowed the closest laser probe volume 
positioning to the blade surface while giving an adequate data rate. The pitch and yaw 
angles which were used for the boundary layer surveys can be found in the summary in 
Appendix C. The boundary layer surveys were performed along a line normal to the blade 
surface at that station. 

d. Wake Surveys 

Wake surveys of various detail were performed at Stations 11, 12, and 13. 
Station 11, the near wake, and Station 13, the far wake, are presented in the next section 
of this report. The laser pitch angle was 5 degrees down for the Station 11 wake surveys. 
The laser was horizontal for the Station 13 surveys. The incremental distance between 
measurement positions was reduced over a series of three surveys at each of the two 
stations, from a foil passage width and large increments to a smaller passage width and 
smaller increments. Again, the laser was aligned prior to each survey using the alignment 
tool. 
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D. COMPUTATIONAL FLUID DYNAMICS 


A quasi-three-dimensional computational fluid dynamics program was used to 
analyze the flow through the cascade. The FORTRAN code. Rotor Viscous Code Quasi- 
3-D (RVCQ3D) Version 300, was written by Rodrick V. Chima of NASA LeRC [Ref 
12]. A 340 x 49 C-type grid was generated using the FORTRAN program “GRAPE” 
[Ref. 13] based on blade manufacturing dimensions [Table 1], Appendix E contains the 
input data used for “GRAPE”. The computational grid generated by “GRAPE” is shown 
in Figure 12. RVCQ3D had the option to use three different turbulence models: the 
Baldwin-Lomax model, the Cebeci-Smith model, and the Wilcox k-omega model; all were 
used. The design inlet-flow angle, Mach number, Reynolds number and pressure ratio 
were other inputs required to run the program. The code was run for various pressure 
ratios until the design inlet flow angle of 36.3 degrees and the Mach number of 0.22 were 
met. [See Appendix E for a sample of the input data used for RVCQ3D], A FORTRAN 
program modified by Williams [Ref. 11] called “PCP” was used to extract blade surface 
pressure distribution data from the solution file in order to plot coefficients of pressure, 
Cp. The program “PCP” is also included in Appendix E. 
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IV. RESULTS AND DISCUSSION 


A. REFERENCE CONDITIONS 

The results of the experimental data will be discussed in the following sections. 
Nominally, all surveys were performed at a Mach number of 0.22 and a Reynolds number 
based on blade chord of 640,000. 


B. BLADE SURFACE PRESSURE DISTRIBUTION 

Figure 13 shows the results of the blade surface pressure distribution 
measurements in terms of the coefficient of pressure, Cp, plotted along the blade chord at 
various positions given by the ratio £/c. Blades 2 and 8 were partially-instrumented with 
eight pressure taps each. The subscript “m” in the figure’s legend denotes the manometer 
readings used on blades 2 and 8. Blade 6 contained 42 pressure taps, including a leading- 
edge tap and a trailing-edge tap and 20 taps on each of the pressure and suction surfaces. 
Two taps were unusable on the pressure surface of blade 6, so a total of 40 pressure 
measurements was made. Gelder’s design intent for Stator 67B [Ref. 1] is shown plotted 
with the solid black line in Figure 13. 

The results from the different blades and different data acquisition methods 
showed good agreement. The Cp calculations from both the water manometer readings 
and the data acquisition system produced very similar results. The Cp results for the blade 
pressure surface were very tightly grouped together for all three blades and both 
measurement techniques. The blade suction-side results showed more loosely grouped 
data for blade numbers 2 and 8 for the two data acquisition methods. Blade number 2 had 
an unusable pressure tap at 0.95 £,/c, and also showed a partially clogged pressure tap at 
0.69 £,/c. The suction side data from blades 2 and 8 corresponded well to the data on 
blade 6 at all other locations. 
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The leading-edge pressure tap on blade number 6 showed the stagnation point to 
be very near the leading-edge, since the leading-edge value was Cp = 0.92. Stagnation 
would give unity. This was a confirmation of the design inlet-flow angle which resulted in 
zero incidence at the leading-edge. The trailing-edge pressure tap showed the value of Cp 
to be 0.10, or very close to zero. 

Gelder’s experimental results for the Stage 67B [Ref. 1] were similar to the 
present results in terms of the relation of the experimental data to the design intent. On 
the pressure side of the blade, Gelder’s experimental Cp data were lower in magnitude 
than the design value, similar to the NPS results. On the suction side of the blade, the 
Gelder experimental data had smaller (negative) magnitudes than design until 
approximately 0.2 £/c, at which point the experimental data had larger magnitudes than 
the design data. 

The experimental Cp profile seemed to show that the design goal of controlled 
diffusion was achieved. On the suction side of the blade, Cp decreased (negatively) very 
sharply from the leading-edge of the blade as the flow accelerated. Cp then decreased 
more gradually until at nearly midchord, where Cp began a near-linear increase, to arrive 
close to zero at the trailing-edge. The near-linear increase showed the diffusion or 
deceleration to be gradual, as the design goal required. 

During tunnel operation, the water manometers connected to the instrumented 
blades showed fluctuations in the pressures. The suction side of the blade had the largest 
fluctuations. The occurrence of fluctuations raised the question of flow separation on the 
suction side. Additional measurements were needed to determine whether the flow had 
separated. 
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c. 


FIVE-HOLE PRESSURE PROBE MEASUREMENTS 


1. Manual Loss and AVDR Calculations 

The manual system was used first to make a downstream survey of 154 mm width, 
with an interval of 2.54 mm. A similar upstream survey was performed. The loss and 
AVDR were calculated using the equations discussed in the previous chapter and found in 
Appendix A. The AVDR was 1.023 and the loss coefficient was 0.029 using the Prandtl 
probe as a reference. The calculated loss coefficient was similar to Gelder’s experimental 
result of 0.030 [Ref 1], 

Figure 14 shows the results of the survey. The static pressure coefficient showed a 
rise downstream of the test section of approximately 1%. The blade wake was well-mixed 
out at the downstream station. The upstream flow angle was approximately 35.5 degrees 
and the downstream flow angle was approximately 2.5 degrees. The nondimensional 
velocity showed a decrease from approximately 0.10 to 0.08. These experimental results 
are slightly different from the LDV results since the measuring station was located 
approximately 2 blade chords upstream and downstream. 

2. HP Automated Data Acquisition System 

The HP automated data acquisition system results which were calculated by 
Armstrong’s “LOSS” program were inconclusive. The calculated loss values ranged from 
-0.013 to 0.011. The AVDR values ranged from 1.04 to 1.05. 

During the surveys, the digital voltmeter showed drift in both the zero and span 
calibration readings. Additionally, data sampling was observed to occur prior to the 
pressure reading from the Scanivalve transducer stabilizing after a step to the following 
channel. The conclusion, following these observations and from the resulting negative 
losses, was that the HP system results were invalid. 
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D. 


LDV RESULTS 


1. Inlet Surveys 

LDV measurements upstream of the test section were performed at Stations 1, 2 
and 3. Stations 1 and 3 will be examined here in order to characterize test section inlet 
flow in both the near and far field. Station 1 was located upstream of the test section at 
30% axial chord (0.30c ac ). Station 1 was surveyed over 167% of passage width, or 254 
mm. Three thousand data points were taken at each position of the survey, with a total of 
41 positions spaced 6.35 mm apart. Results at Station 1 in the form of velocity ratios 
referenced to the inlet velocity condition, V re f, turbulence percentage referenced to V re f, 
and the Reynolds stress correlation coefficient, Cuv, are plotted in Figure 15. 

Station 1 total velocity ratio, W/V re f, was nearly uniform across the passage span. 
Both the axial velocity ratio, U/V re f, and the tangential velocity ratio, V/V re r, showed a 
slight variation across the passage. The potential influence of the blades was felt as far 
upstream as 30% axial chord, which resulted in the depressions in velocity spaced one 
blade passage width apart. The axial turbulence, T u , was measured to be 
uncharacteristically high for the position in the flow, while the tangential turbulence, T v , 
was consistent at 2%. The most likely reason for the high turbulence of T u was due to the 
laser optics misalignment, which was resolved during later surveys. The Reynolds-stress 
correlation coefficient was below 0.1, showing the flow to be random or uncorrelated. 

The survey was repeated at Station 1 over a 165.1 mm width at 27 positions with 
6.35 mm interval spacing, or 108% of passage width due to the turbulence discrepancy 
noted above, with the results shown in Figure 16. The mean velocities showed as much 
variation as was recorded previously. Both turbulence percentages were now of 
equivalent magnitude, indicating that the optics problem from the first survey was 
resolved. The correlation coefficient was again less than 0.1. 

Station 3 was located upstream 5% of an axial chord (0.05c ac ) from the leading- 
edge. Two surveys were performed. For the first survey with the laser horizontal, the 
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survey width was only 91% of blade passage width or 142 mm, and 72 data points spaced 
2 mm apart were taken. The survey did not cover the full passage due to blade 
interference with the two laser beams. The survey was repeated at Station 3 over a width 
of 203.2 mm, with 41 positions at a spacing of 5.08 mm, with the laser pitched upward 5 
degrees. Figure 17 shows the overlay of the data of the two surveys, with and without 
laser pitch. The figure shows the repeatability of the velocity, turbulence and correlation- 
coefficient data for the two surveys was excellent. The proximity of the leading edge of 
the blades is evident in all three plots, shown by the velocity gradients near y/S = 0 and 1. 
The total velocity decreased to a minimum along the stream line leading to the stagnation 
point on the leading-edge of the blade, and increased dramatically on either side as the 
flow proceeded around the leading-edge. The maximum total velocity was to the right of 
the stagnation point, indicating the suction side of the blade. Periodicity was also shown 
by the repetition of the velocity distributions. Turbulence remained constant at 
approximately 2%, with slight increases near the stagnation streamlines. The correlation 
coefficient showed a peak near the stagnation points of the leading-edge on the pressure 
side of the blades, and remained less than 0.2 throughout the survey, indicating a slight 
reorientation of the turbulence around the leading-edge. 

2. Passage Surveys 

Passage surveys were conducted at Stations 4 through 10 in the passage between 
blades 3 and 4 with neither laser yaw nor laser pitch used during these surveys. Flow at 
Stations 5, 7, 8, 9 and 10 will be discussed in this section, and all reduced passage data 
can be found in Appendix C. 

Figure 18 shows the Station 5 results. The velocity ratios were smooth and 
showed a decrease across the passage from the suction side to the pressure side of the 
passage. The turbulence percentages were nearly equal and maintained a value of 
approximately 2%. The correlation coefficient ranged from approximately -0.1 to 0.2. 
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The Station 7 survey results are presented in Figure 19. The velocity data in the 
figure showed a decrease for two points to the left of the maximum in both total and axial 
velocity on the suction side of the blade. The velocities then decreased gradually as the 
distance away from the suction surface increased. Turbulence peaked at 12% in the axial 
direction and decreased to 2%. Turbulence remained at 2% in the tangential direction. 
The correlation coefficient ranged from 0 to 0.2. These results indicated that the passage 
survey for Station 7 reached slightly into the boundary layer, which suggested that the 
boundary layer was growing in thickness, compared with Station 5. 

Station 8 results are presented in Figure 20. The velocity ratios were smooth, with 
a slight decrease in magnitude as y/S increased. Turbulence was approximately 2 %, 
similar to previously discussed stations. The correlation coefficient ranged from 0.1 to 
0.3. 

Figure 21 shows the results from the survey at Station 9. The velocity ratios 
remained smooth, turbulence remained at 2%, and the correlation coefficient again ranged 
from approximately 0.1 to 0.3. Figure 22 shows Station 10 results. The velocity ratios 
were similar to those at Station 9, the turbulence remained approximately 2%, and the 
correlation coefficient ranged from 0.1 to 0.3, as before. 

3. Wake Surveys 

Wake surveys were performed at Stations 11, 12 and 13. The results from 
Stations 11 and 13, the near and far wake, will be discussed. A series of three surveys 
were performed at each station, ranging in extent from a full passage width and decreasing 
to the width of just the wake. 

Figure 23 shows Station 11 results for nearly the full passage width. The laser 
remained horizontal for this survey, so blade interference with the laser beams prevented a 
full passage-width survey. The velocity ratios are shown to be nearly uniform. 
Turbulence remained approximately 2% until the wake was reached, at which point the 
turbulence began to increase. The correlation coefficient remained below 0.3. 
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The laser was pitched 5 degrees down for the next two wake surveys at Station 
11. Figures 24 and 25 show the results. The total velocity ratio decreased and then 
increased again as the wake was traversed during the survey. The axial velocity ratio was 
similar to the total velocity ratio. The tangential velocity ratio showed an initial increase 
and then decreased below zero as the wake was surveyed. The wake turbulence showed 
two distinct peaks in the axial direction as the wake was traversed, while the tangential 
direction showed a single peak. The maximum axial turbulence was 18%, and the 
maximum tangential turbulence was 16%. The correlation coefficient started at a 
magnitude of 0.2 and became -0.2 as the wake was traversed, and returned to a value of 
approximately 0. Figure 26 shows the results of the two Station 11 wake surveys plotted 
together. The data of the two surveys correlated very well, as shown by the close overlay. 

Station 13 was surveyed over three different ranges from an entire blade passage 
width to just the width of the wake. Figure 27 shows the first survey over more than the 
passage width. The velocity ratio showed a decrease through the wake, giving an 
indication of the necessary width for the following survey. The turbulence was 
approximately 2% over the passage, until the wake was encountered. In the wake, 
turbulence increased to approximately 14%. The correlation coefficient decreased from 
0.2 to -0.2 through the wake, and remained approximately 0.1 outside of the wake. 

Figures 28 and 29 show more detailed wake surveys at Station 13. The velocity 
ratios again showed a decrease as the wake was traversed, and then an increase. The axial 
turbulence in the wake showed two peaks, similar to the Station 11 results, and had a 
maximum magnitude of 14%. The tangential turbulence showed a single maximum peak, 
similar to the Station 11 results, with a maximum magnitude of 12%. The correlation 
coefficient went from 0.2 to -0.2 in the wake, while remaining at approximately 0.1 
outside of the wake. Figure 30 shows the two Station 13 wake surveys plotted together. 
The data correlated well and only the axial turbulence, T u> showed a slight difference in 
value on the right side of the wake. 
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The axial turbulence results in Figures 26 and 30 are of interest. While both 
stations show two peaks in T u , Station 11 showed the maximum peak to be located on the 
left side of the wake, while Station 13 showed the maximum magnitude to be located on 
the right hand side. While asymmetry in the wake turbulence was measured in the wake of 
the Stator 67A cascade [Ref. 5], a reversal of the asymmetry from one side to another was 
not previously observed. ' 

4. Boundary Layer Surveys 

Boundary layer surveys were performed at four stations, normal to the surface of 
the blade at the specified station. Station 5, 7 and 9 were surveyed on the suction side, 
while Station 8 was surveyed on the pressure side. A combination of laser pitch and yaw 
was used in order to position the LDV probe volume as close to the blade as possible. 

A total of three surveys were performed at Station 5 with three different 
combinations of laser yaw and pitch. The combination which worked best was a yaw of 5 
degrees and a pitch of 0 degrees, due to beam reflection from the cascade wall and off the 
brass shims on the far endwall. The results of the boundary-layer survey at Station 5 are 
plotted in terms of the ratio of distance from the wall and blade chord, and are shown in 
Figure 31. The velocity ratios showed only a slight decrease on the left side, at the 
minimum y/S value, which corresponded to the edge of the boundary layer. The 
turbulence was a maximum for approximately 4 survey points at the boundary layer edge, 
which decreased to approximately 2% away from the boundary layer. The correlation 
coefficient was a maximum near the blade surface, and decreased to 0.1. No data could be 
taken close to the blade. 

The results of the boundary-layer survey taken at Station 7 are shown in Figure 32. 
The velocity ratios showed a low flat area near the blade, which then rose to near 
freestream values. The axial turbulence started near 12%, rose to a peak of 40% and then 
began to decrease near the freestream. Tangential turbulence remained approximately 2%. 
The correlation coefficient fluctuated from -0.1 to 0.1 for the survey. The high level of 
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axial turbulence suggested that the flow had separated in the vicinity of Station 7. The 
boundary layer data, along with the fluctuating pressure noted during Cp measurements, 
and the passage survey results, collectively indicated flow separation. 

The results from the Station 8 pressure-side boundary-layer survey are shown in 
Figure 33. The abscissa was plotted in reverse order, since the pressure side of blade 4 
was on the right-hand side of the survey passage. The velocity ratios showed a decrease 
in the total velocity for the last seven points of the survey, on the right side. The 
turbulence increased for the seven positions near the blade, to a maximum of 6%. The 
correlation coefficient showed an increase from approximately 0.3 to 0.4 in the boundary 
layer. 

Station 9 boundary layer results are plotted in Figure 34. This survey was 
performed on the suction side of the blade. The velocity ratio showed a decrease for the 
ten points near the blade surface, and the turbulence increased in both the axial and 
tangential direction as the blade surface was approached. The boundary layer profile 
suggested that the flow had reattached itself to the suction surface prior to Station 9. The 
correlation coefficient decreased below 0 in the boundary layer, and then increased to 0.1 
in the passage. 


5. Combination of Boundary Layer and Passage Surveys 

Two Station 7 survey results are presented together in Figure 35. The passage 
survey at Station 7 indicated the boundary layer was thicker than at other stations, 
allowing a few data points to be taken inside the boundary layer without laser pitch or 
yaw. The boundary layer survey at Station 7, performed along a path normal to the blade 
surface at Station 7, showed that a large region of separation existed. In order to study 
the flow close to the blade surface with the same reference as the passage surveys, the 
laser was yawed to perform another passage survey at Station 7. Data from this survey 
with the laser yawed was plotted with the Station 7 passage survey data [from Figure 19], 
in Figure 35. The velocity data in the figure show a flat region for both total and axial 
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velocity in the boundary layer, and then a rapid rise to freestream conditions. Turbulence 
peaked at 24% in the axial direction, and remained at 2% in the tangential direction. The 
correlation coefficient was erratic in the boundary layer, ranging from -0.1 to 0.1. The 
correlation coefficient then settled to a positive value of 0.2 outside the boundary layer. 

Hobson and Shreeve [Ref. 9] located separation near the leading edge of Stator 
67A blades in the NPS cascade. The velocity profile at Station 3 from Hobson and 
Shreeve [Ref. 9] had a similar shape to the Stator 67B Station 7 velocity profile [Figure 
35], The turbulence results from their measurements and those at Station 7 are quite 
similar. The results of the combined data in Figure 35, Cp observations, and the boundary 
layer data all indicate that flow separation was present on the suction side at Station 7. 

The passage survey data and the boundary layer data at Stations 8 and 9 are 
plotted in Figures 36 and 37, respectively. The results were of interest since the passage 
survey was nearly normal to the blade at those stations. The boundary layer extent was 
evident in these figures. The velocity ratios decreased in the boundary layer and the 
turbulence increased significantly. 


E. CFD RESULTS 

Rotor Viscous Code Quasi-3-D was run using three turbulence models. The 
pressure ratio, prat, was iterated upon until the inlet-flow angle converged to within 0.02 
degrees. The CFD predictions are plotted in Figure 38 together with the experimental 
results for the pressure distribution. Table 3 gives the values of the pressure ratio, prat, 
which yielded design inlet-flow angle when the CFD code was used with the three 
different turbulence models. 

All computed results showed a sharp acceleration on the suction side initially, and 
then a gradual deceleration for most of the rest of the suction side. Near £/c = 0.70 on the 
suction side, the code predicted a slight re-acceleration. These CFD results for the suction 
side did not match the experimental data which showed a more gradual acceleration. 
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peaking at approximately 5,/c = 0.40, and then a near-linear deceleration. The predicted 
results on the pressure side of the blade compared quite well with the experimental data. 

The Wilcox k-omega turbulence model had the capability to modify the freestream 
turbulence, which was computationally changed from 0.02 to 0.04. The result from an 
increase in turbulence from 0.02 to 0.04 showed no significant change in the predicted Cp 
results. 
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Figure 14. Manual Loss Survey Results. 
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Figure 16. Repeat Station 1 Survey Results for 108% Passage Width. 
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Figure 17. Station 3 Survey Results With and Without Laser Pitch. 
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Figure 19. Station 7 Passage Survey Results. 
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Figure 20. Station 8 Passage Survey Results. 
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Figure 22. Station 10 Passage Survey Results. 
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Figure 23. Station 11 Survey Results. 


50 








Ratio 



Station 11 Wake Velocity Ratios 



15 0.20 0.25 0.30 0.35 


y/S 


Station 11 Wake Correlation Coefficient 



Figure 24. Station 11 Wake Survey Number 1 Results. 
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Figure 25. Station 11 Wake Survey Number 2 Results. 
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Figure 26. Station 11 Wake Surveys 1 and 2 Results. 
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Figure 27. Station 13 Survey Results. 
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Figure 29. Station 13 Wake Survey Number 2 Results. 
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Figure 30. Station 13 Wake Surveys 1 and 2 Results 















Figure 31. Station 5 Boundary Layer Survey Results. 
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Figure 33. Station 8 Boundary Layer Survey Results. 
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Figure 35. Station 7 Boundary Layer and Passage Survey Results. 
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Figure 36. Station 8 Boundary Layer and Passage Survey Results. 
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Model 

Cebeci-Smith 

Baldwin-Lomax 

Wilcox k-omega 

Pressure Ratio (prat) 

0.9765 

0.97635 

0.97625 


Table 3. Pressure Ratios for Turbulence Models. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Compressor Stator 67B cascade blading was successfully installed and 
experimentally tested at design conditions in the NPS Low-Speed Cascade Wind Tunnel. 
The experiments were conducted at a design inlet-flow angle of 36.3 degrees, a Mach 
number of 0.22 and a Reynolds number based on chord of 640,000. 

Experimental blade surface pressure measurements were obtained using water 
manometers and an HP automated data acquisition system. The resulting Cp plots 
showed that the flow on the suction side rapidly accelerated from the leading edge, more 
gradually accelerated until midspan, and then gradually decelerated to the trailing edge. 
Fluctuations in pressure on the suction side of the blades observed during the experiments 
suggested unsteady flow or separation had occurred on the blade, and this was later 
confirmed with LDV passage and boundary layer surveys. 

Experimental measurements yielded a loss of 0.029. These results compared well 
to NASA LeRC loss values of 0.030 [Ref. 1], However, the NASA results were from 
tests of the entire stage, so the comparison is not totally valid. 

The LDV results characterized the flow in the inlet, in the blade passage, in the 
boundary layers and at the exit of the test section at design conditions. Inlet surveys 
measured the upstream influence of the blades. Passage surveys showed a smooth 
acceleration up to a separation point in the vicinity of Station 7 on the suction side of the 
blade, and then a smooth deceleration and reattachment of the flow. With the presence of 
separation, although the blade produced a useful force, the controlled diffusion design 
goal of preventing separation at design incidence was not achieved. Boundary-layer 
surveys were successfully performed at three suction-side stations, and one pressure-side 
station. The boundary-layer surveys defined the growth of the boundary layer throughout 
the passage. Wake surveys showed the influence of mixing in the wake. 
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CFD predictions of the blade surface pressure coefficients using RVCQ3D Version 
300 did not correlate well with the experimental results. Three different turbulence 
models were used, with similar results. 


B. RECOMMENDATIONS 

More extensive boundary-layer surveys and analysis should be performed at design 
incidence to determine more exactly the regions of separation and reverse flow. Laser 
rotation, in addition to laser pitch and yaw, could allow the LDV probe volume to be 
positioned closer to the blade surface at the measurement stations. Flow visualization 
using a laser sheet should be performed to help determine the separation location, and to 
determine the extent of the separation region. 

The five-hole probe data acquisition system should be improved by updating the 
survey equipment, and controlling software. A computer-controlled automatic-traverse 
mechanism for the five-hole probe, and an automated yaw system to align the five-hole 
probe with the flow, would permit surveys to be accomplished more rapidly. An update of 
the HP BASIC software to include self-adjusting time pauses to allow settling of pressure 
data prior to recording, point-by-point data recording in case of system failure, and an 
option to use Prandtl probe pressures for reference should be implemented immediately. 
The option of using software based on “LabVIEW” to control the system should be 
investigated. Finally, the five-hole probe calibration coefficients should be verified before 
further probe surveys are conducted. 

Additional analysis with RCVQ3D should be performed to try to reproduce the 
experimental results. Different grid sizes and number of iterations should be investigated. 
Also, the incidence angle should be varied in an attempt to match the inlet flow conditions. 
Loss calculations should be included in the CFD analysis. Experimentation with a full 
three dimensional code might provide better results. 
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When the flow at the design incidence angle is initially characterized, positive off- 
design incidence angles should be set to determine immediately the available range of 
incidence and the on-set of stall. This will allow the future test plan to be determined most 
optimally. 
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APPENDIX A. FIVE-HOLE PROBE EQUATIONS 


Dimensionless velocity coefficient: p = ——— 

P. 

Average yaw pressure: p 23 = — * 

Pitch angle coefficient: y^, = ——— 

P. “ P23 


Dimensionless velocity coefficients [Ref. 8] Cij: 



Cil 

C i2 

c i3 

Ci4 

c i5 

C.j 

0.015926 

4.932133 

-153.66876 

3137.9614 

-24299.005 

c 2j 

-0.003563 

0.699505 

-62.977261 

2068.3721 

-20872.148 

Csj 

0.080098 

-24.844173 

1980.4954 

1980.4954 

541835.5 


Dimensionless velocity polynomial: 

X = {c„ +c, ; p+c„ p ! +C ]4 p 5 +c„p 4 } 

+ {c 2 l + C 22 P + C 23 P +C 24 "P +C 25 -P |■ y 5_hole 

+ {c 3 l + C32 P + C J3 P +CJ4-P +C 3 5'P | ' y 5-hole 


Reference dimensionless velocity: 



Reference flow function: 


Axial-Velocity-Density Ratio: 


K = 


‘‘ref L 


I - X 2 


1— X; 


ref 


y-i 


CO; 


5-hole) 


AVDR = 


f K - d * 

f K ^ 
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Dimensionless total pressure: 


Dimensionless static pressure: 


Loss coefficient: 


P, 

P, 


p _ ^23 

W* p 


CO = 


f C Dtus K us dx--— f C Dtds K ds dx 

Jo ptus us AVDR **° p 


fCpu^dx - j o C psus K us dx 
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APPENDIX C. LDV SUMMARY AND REDUCED DATA 


|LDV Summary 
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| Station 7- Repeat Boundary Layer Survey 
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APPENDIX D. REFERENCE VELOCITY CODE 


FORTRAN CODE “CALIBI” 


C 


C 

PROGRAM TO COMPUTE THE CALIBRATION CURVE FOR THE NPS LOW SPEED C 
CASCADE WIND TUNNEL. 

C 

c 

- c 


A STRAIGHT LINE IS FITTED THROUGH THE REFERENCE CONDITIONS OF 
THE TUNNEL AT DIFFERENT SPEEDS. 

THE REFERENCE VELOCITY IS THEN OBTAINED, BY NEWTON'S METHOD, 
DEPENDING ON THE TUNNEL PLENUM PRESSURE AND TEMPERATURE 


- c 

PROGRAM CALIBRATE 
IMPLICIT REAL*8(A-H,0-Z) 

PARAMETER (NP=6) 

DIMENSION VA(NP),VT(NP),PA(NP),PP(NP),TP(NP) 

DIMENSION VTOT(NP),X(NP),ANUX(NP),PR(NP) 

DIMENSION PAR( 100),PPR( 100),TPR( 100),PRR( 100), ANUXR( 100) 

DIMENSION VREF(IOO) 

CHARACTER* 14 NAME(IOO) 

OPEN(UNIT=10,FILE='CALIB.DAT',STATUS='UNKNOWN’) 

OPEN(UNIT= 11 ,FILE='REFER.DAT',STATUS='UNKNOWN') 
OPEN(UNIT=12,FILE='CALIB.OUT',STATUS-UNKNOWN') 

PRINT BANNER 

WRITE(*,500) 

WRITE( 12,500) 

DO FORMAT(1 X,78('C'),/1 X,’C',76X,'C7, 

#1X,’C',20X,'OUTPUT FROM PROGRAM CALIBRATE',27X,'C'/1X,'C',76X,'C'/ 
#1X,78('C')// 

#5X,'LEAST SQUARES STRAIGHT LINE CURVE FIT IS USED'/ 

#5X,TO DETERMINE TUNNEL CHARACTERISTICS AT DIFFERENT SPEEDS'// 

#5X,'NEWTON S METHOD IS USED TO DETERMINE THE REFERENCE VELOCITY'/ 
#5X,'FROM THE RECORDED AMBIENT PRESSURE AND TUNNEL PLENUM'/ 
#5X,'PRESSURE AND TEMPERATURE'/) 

INITIALIZE AIR, MERCURY AND WATER PROPERTIES 
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c 

RHOW=1000.D0 
RHOHG=13000.D0 
RHOA= 1.2256D0 
CPA=1005.D0 
GAMMA=1.4D0 
GM1=GAMMA-1.D0 
C 

C BEGIN DETERMINING TUNNEL CHARACTERISTICS 
C 

WRITE(*,501) 

WRITE(12,501) 

501 FORMAT(//1X,'BEGIN DETERMINING TUNNEL CHARACTERISTICS’,/ 

# IX,’FROM THE FOLLOWING MEASURED VALUES’,// 

#3X,'AXIAL VEL. TANGENTIAL VEL. AMBIENT PRESS.' 

#3X,'PLENUM PRESS. PLENUM TEMP.'/ 

#3X,'M PER SEC. M PER SEC. INCHES MERCURY’ 

#3X,'INCHES WATER DEG. C.7) 

C 

C READ IN DATA POINTS 
C 

C VA = AXIAL VELOCITY - MEASURED BY LDV (M/S) 

C VT = TANGENTIAL VELOCITY - MEASURED BY LDV (M/S) 

C PA = AMBIENT PRESSURE (INCHES MERCURY) 

C PP = PLENUM PRESSURE (INCHES WATER) 

C TP = PLENUM STAGNATION TEMPERATURE (DEG. C.) 

C 

DO 1 1=1,NP 

RE AD( 10,100) VA(I), VT (I),PA(I),PP(I),TP(I) 

WRITEC", 110) VA(I), VT(I),P A(I),PP(I),TP(I) 

WRITE(12,110)VA(I), VT(I),PA(I),PP(I),TP(I) 

1 CONTINUE 

100 F0RMAT(1X,5F8.4) 

110 FORMAT(3X,F8.4,7X,F8.4,8X,F8.4,7X,F8.4,7X,F8.4) 

C 

C CALC VTOT, X, ANUX AND PR 
C 

WRITER,510) 

WRITE(12,510) 

510 FORMAT(/5X,'CALCULATED VALUES FOR THE TUNNEL CONFIGURATION',// 
#3X,TOTAL VELOCITY’,6X, 'MACH NUMBER',9X,'MACH NUMBER FUNCT.',3X, 
#'PRESSURE RATIO'//) 

C 

DO 2 1=1,NP 

VTOT(I)=DSQRT(VA(I)**2+VT(I)**2) 

X(I)=VTOT(I)/DSQRT(2.DO*CPA*(TP(I)+273.16DO)) 

ANUX(I)=(G AMMA/GM 1)*(X(I)**2)*(1.D0-(X(I)*»2))**(1 .D0/(GM 1)) 
PR(I)=1.DO-(RHOHG*PA(I))/(RHOW*PP(I)+RHOHG*PA(I)) 

WRITEC", 101) VTOT(I),X(I),ANUX(I),PR(I) 

WRITE( 12,101) VTOT (I),X(I), ANUX(I),PR(I) 

2 CONTINUE 
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101 FORMAT(4D20.11) 

C 

C CALL THE LEAST SQUARES SUBROUTINE TO FIT A STRAIGHT LINE THROUGH THE 
C DATA; X-AXIS = PR, Y-AXIS = ANUX (MACH NO. PARAMETER) 

C 

WRITE(*,520) 

WRITE(12,520) 

520 FORMAT(/10X,'CALLING LEAST SQUARES SUBROUTINE',/IX, 

#'TO DETERMINE THE PRESSURE RATIO AS A FUNCTION OF MACH NO. PARAM' 
mx; 'PRESSURE RATIO = A1 * ANUX + AO',//) 

c 

CALL LEASTSQUARE(NP,ANUX,PR,A0,A1) 

C 

WRITE(*,530)Al,AO 

WRITE(12,530)A1,A0 

530 FORMAT(/10X,'A1 = ',D20.11,’ AO = ',D20.11/) 

C 

C READ IN REFERENCE CONDITIONS 
C 

WRITE(*,199) 

WRITE(12,199) 

199 F0RMAT(/1X,'REFERENCE CONDITIONS FOR EACH RUN'/ 

# IX,'AMBIENT PRESSURE PLENUM PRESSURE PLENUM TEMPERATURE', 

#3X,'RUN NAME'/ 

mx; INCHES MERCURY INCHES WATER DEGREES CELSIUS'//) 

C 

C THE NUMBER OF EXPERIMENTS 
C 

READ(11,198)NE 
198 F0RMAT(1X,I3) 

C 

C THE REFERENCE CONDITIONS 
C 

DO 99 1=1,NE 

READ(11,200)PAR(I),PPR(I),TPR(I),NAME(I) 

WRITE(*,201)PAR(I),PPR(I),TPR(I),NAME(I) 

WRITE( 12,201 )P AR(I),PPR(I),TPR(I),N AME(I) 

99 CONTINUE 

200 FORMAT( 1X,3F8.4,A 14) 

201 FORMAT(2(6X,F8.4),10X,F8.4,11X,A14) 

C 

C CALCULATE THE PRESSURE RATIO 
C 

DO 88 1=1,NE 

PRR(I)= 1 .DO-(RHOHG*PAR(I))/(RHOW*PPR(I)+RHOHG*PAR(I)) 

C 

ANUXR(I)=(-A0+PRR(I))/A 1 
C 

WRITE(*,550)I,PRR(I),ANUXR(I),NAME(I) 

WRITE( 12,550)1, PRR(I), ANUXR(I),NAME(I) 

550 FORMAT(/lX,'I = ',13,/ 
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#' PRESSURE RATIO = ',F10.5,' MACH NUMBER PARAMETER = \E12.4,/ 

#' RUN NAME = ',A14/) 

C 

C NEWTON METHOD TO DETERMINE THE ROOTS OF THE EQUATION FOR ANUX 
C 

C GUESS INITIAL VALUE OF VREF = VTOT(5) 

C 

WRITE(*,560) 

WRITE( 12,560) 

560 FORMAT(//10X,'BEGIN NEWTON ITERATION'/) 

C 

KOUNT=l 

VREF(I)=VTOT(5) 

XR=VREF(I)/DSQRT(2.D0*CPA*(TPR(I)+273.16D0)) 
1000F=/GAMMA/GM1)*(XR**2)*/1.D0-(XR**2))**(1.D0//GM1))-ANUXR(I) 
DFDX=/GAMMA/GM1)*2.D0*XR*((1.D0)-(XR**2)*(1.D0/(GM1)))* 

# (1 .D0-(XR**2)*(1 .DO/GM1 )*(1 ,D0/(1 ,D0-(XR**2)))) 

TERM=F/DFDX 

C 

WRITE/*,570)KOUNT,XR, TERM 
WRITE/12,570)KOUNT,XR,TERM 

570 FORMAT(1 X,'ITERATION NUMBER ',12,' MACH NO. PARAM. = ',F8.6, 

#' ERROR TERM = \D12.4) 

C 

XR=XR-/F/DFDX) 

KOUNT=KOUNT +1 
IF(DABS/TERM).LT.1.D-10)GO TO 999 
IF/KOUNT.LT. 10.OR.DABS/TERM).GT. 1 D-10)GO TO 1000 
C 

C CALCULATE THE REFERENCE VELOCITY 
C 

999 VREF/I)=XR*DSQRT(2.D0*CPA*(TPR(I)+273.16)) 

WRITE(*,580)VREF(I) 

WRITE/12,580) VREF/I) 

580 FORMAT(/1X,’VREF = \F20.11/) 

88 CONTINUE 
C 

C PRINT FINAL RESULT 
C 

WRITE/12,585) 

585 FORMAT(/lX,'EXPERIMENT NUMBER REFERENCE VELOCITY NAME'/) 

DO 777 I=1,NE 

WRITE/12,590)1, VREF/I),NAME/I) 

777 CONTINUE 

590 FORMAT(8X,I2,16X,F8.4,10X,A14) 

STOP 

END 

c-c 

c 
c 

C A LEAST SQUARES CURVE FIT FOR A STRAIGHT LINE THROUGH NOISY 
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C DATA 
C 

C ALGORITHM TAKEN FROM: NUMERICAL METHODS, ROBERT W. HORNBECK 
C PAGES 122 - 130, QUANTUM, 1975 

C 


SUBROUTINE LEASTSQUARE(N,X,FX,A0,A1) 

IMPLICIT REAL* 8(A-H,0-Z) 

DIMENSION COEFF(2,2),RHS(2),X(N),FX(N) 

C 

C SET UP COEFFICIENT MATRIX AND RIGHT HAND SIDE. 
C 

COEFF(I,l)=N 

COEFF(1,2)=O.DO 

COEFF(2,2)=O.DO 

C 

RHS(1)=0.D0 

RHS(2)=0.D0 

C 

DO 2 1=1,N 

COEFF( 1,2)=COEFF( 1,2)+X(I) 
COEFF(2,2)=COEFF(2,2)+X(I)**2 
RHS( 1 )=RHS( 1 )+FX(I) 

RHS(2)=RHS(2)+X(I)*FX(I) 

2 CONTINUE 
COEFF(2,1 )=COEFF( 1,2) 

C 

C PRINT MATRIX EQUATION 
C 

WRITE(*,110) 

WRITE(12,110) 

110 FORMAT(/5X,'MATRIX EQUATION',/) 

DO 3 1=1,2 

WRITE(*,102)(COEFF(I,J),J=1,2),I-1,RHS(I) 

WRITE(12,102)(COEFF(I,J),J=1,2),I-1,RHS(I) 

3 CONTINUE 

102 FORMAT(5X,2(2X,D 10.2),' A',I1,3X,D10.2) 

C 

C GAUSS ELIMINATION 
C 

TERM=COEFF(2, l)/COEFF( 1,1) 

C 

COEFF(2,2)=COEFF(2,2)-COEFF( 1,2)*TERM 
RHS(2)=RHS(2)-RHS(1)*TERM 
C 

Al=RHS(2)/COEFF(2,2) 

A0=(RHS( 1 )-COEFF( 1,2)* A1 )/COEFF( 1,1) 

C 

RETURN 

ENDD 
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FORTRAN INPUT FILE “REFER.DAT 


28 

29.9499 12.1000 21.1111 0502sl3a. 

29.9499 11.9000 21.1111 0502sl3b 

29.9499 11.9000 21.1111 0502sl3c 

29.9499 11.9000 21.1111 0502sl2a 

29.9499 11.9000 21.1111 0502sl2b 
29.8684 12.0000 20.0000 0506slla 

29.8481 12.1000 20.5556 0506sl0_ 

29.8481 12.1000 20.5556 0506s9_ 

29.8481 12.1000 20.5556 0506s8_ 

29.8481 11.9000 21.6667 0506s7a_ 
30.0517 12.0000 21.1111 0509s6a_ 
30.0313 12.1000 20.5556 0509s5_ 

29.8888 12.0000 21.1111 0509s4_ 

29.8888 11.9500 20.5556 0507s3a_ 

29.8888 11.9000 20.0000 0507s2a_ 
29.9702 12.0000 21.1111 0401sla_ 

29.9295 11.9000 21.1111 0520s6b_ 

29.9295 11.8000 22.2222 0520sllb 

29.9295 11.9500 22.2222 0520sllcf 
29.9091 12.2000 22.2222 0520bp8_ 

29.9499 11.9500 20.0000 0526s2b_ 

29.9499 12.0000 20.5556 0526s3b_ 

29.9295 11.9500 20.5556 0603bs5_ 
29.9091 12.0000 25.5556 0708bs9_ 

29.9295 11.9000 23.3333 0404slb_ 
29.8684 12.0000 22.2222 0901bs7_ 

29.8481 11.8500 24.4444 0901bs7b_ 

29.8481 11.8000 25.0000 0901s7b 
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FORTRAN OUTPUT FILE “CALIB.OUT 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c c 

C OUTPUT FROM PROGRAM CALIBRATE C 

C C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

LEAST SQUARES STRAIGHT LINE CURVE FIT IS USED 
TO DETERMINE TUNNEL CHARACTERISTICS AT DIFFERENT SPEEDS 

NEWTON S METHOD IS USED TO DETERMINE THE REFERENCE VELOCITY 
FROM THE RECORDED AMBIENT PRESSURE AND TUNNEL PLENUM 
PRESSURE AND TEMPERATURE 


BEGIN DETERMINING TUNNEL CHARACTERISTICS 
FROM THE FOLLOWING MEASURED VALUES 


AXIAL VEL. TANGENTIAL VEL. AMBIENT PRESS. PLENUM PRESS. PLENUM TEMP. 
M PER SEC. M PER SEC. INCHES MERCURY INCHES WATER DEG. C. 


33.4750 

24.5991 

29.9499 

3.6000 

20.5556 

42.7350 

31.4532 

29.9499 

5.8000 

21.6667 

50.4851 

37.3610 

29.9499 

8.1000 

22.7778 

56.1076 

41.5132 

29.9295 

10.1000 

23.3333 

60.7966 

45.0193 

29.9295 

11.9000 

23.3333 

65.9383 

48.5611 

29.9295 

14.0000 

23.8889 


CALCULATED VALUES FOR THE TUNNEL CONFIGURATION 
TOTAL VALOCITY MACH NUMBER MACH NUMBER FUNCT. PRESSURE RATIO 


0.41541441306E+02 
0.53062076997E+02 
0.62805968212E+02 
0.69795476587E+02 
0.75650273919E+02 
0.81890413603E+02 


0.54065475267E-01 
0.68929132229E-01 
0.81433418389E-01 
0.90411111093E-01 
0.97995252042E-01 
0.10597930366E+00 


0.10I56165233E-01 
0.16432468008E-01 
0.22827031824E-01 
0.28028521818E-01 
0.32809629573E-01 
0.38216120664E-01 


0.9161501234 IE-02 
0.14678018680E-01 
0.20379988663E-01 
0.25301646987E-01 
0.29677031599E-01 
0.34732257708E-01 


MACH NO. PARAM 


CALLING LEAST SQUARES SUBROUTINE 
TO DETERMINE THE PRESSURE RATIO AS A FUNCTION OF 

PRESSURE RATIO = A1 * ANUX + A0 


MATRIX EQUATION 
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0.60E+01 0.15E+00A0 0.13E+00 

0.15E+00 0.42E-02A1 0.38E-02 

Al= 0.91276427032E+00 AO = -0.26460149149E-03 


REFERENCE CONDITIONS FOR EACH RUN 

AMBIENT PRESSURE PLENUM PRESSURE PLENUM TEMPERATURE RUN NAME 
INCHES MERCURY INCHES WATER DEGREES CELSIUS 


29.9499 

12.1000 

21.1111 

0502s13a 

29.9499 

11.9000 

21.1111 

0502sl3b 

29.9499 

11.9000 

21.1111 

0502sl3c 

29.9499 

11.9000 

21.1111 

0502sl2a 

29.9499 

11.9000 

21.1111 

0502sl2b 

29.8684 

12.0000 

20.0000 

0506s11 

29.8481 

12.1000 

20.5556 

0506s10 

29.8481 

12.1000 

20.5556 

0506s9 

29.8481 

12.1000 

20.5556 

0506s8 

29.8481 

11.9000 

21.6667 

0506s7 

30.0517 

12.0000 

21.1111 

0509s6a 

30.0313 

12.1000 

20.5556 

0509s5 

29.8888 

12.0000 

21.1111 

0509s4 

29.8888 

11.9500 

20.5556 

0507s3a 

29.8888 

11.9000 

20.0000 

0507s2a 

29.9702 

12.0000 

21.1111 

0401sla 

29.9295 

11.9000 

21.1111 

0520s6b 

29.9295 

11.8000 

22.2222 

0520slla 

29.9295 

11.9500 

22.2222 

0520sllb 

29.9091 

12.2000 

22.2222 

0520bp8 

29.9499 

11.9500 

20.0000 

0526s2b 

29.9499 

12.0000 

20.5556 

0526s3b 

29.9295 

11.9500 

20.5556 

0603bs5 

29.9091 

12.0000 

25.5556 

0708bs9 

29.9295 

11.9000 

23.3333 

0404slb 

29.8684 

12.0000 

22.2222 

0901bs7 

29.8481 

11.8500 

24.4444 

0901bs7b 

29.8481 

11.8000 

25.0000 

0901s7b 

EXPERIMENT NUMBER 

REFERENCE VELOCITY NAN 

1 

75.9548 

0502s 13a_ 


2 

75.3334 

0502s13b 


3 

75.3334 

0502s13c 


4 

75.3334 

0502s12a_ 


5 

75.3334 

0502s12b 


6 

75.6033 

0506s lla_ 


7 

76.0105 

0506s10 


8 

76.0105 

0506s9 
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9 

76.0105 

0506s8_ 

10 

75.5311 

0506s7a_ 

11 

75.5184 

0509s6a_ 

12 

75.7816 

0509s5_ 

13 

75.7209 

0509s4_ 

14 

75.4939 

0507s3a_ 

15 

75.2668 

0507s2a 

16 

75.6195 

0401sla 

17 

75.3587 

0520s6b 

18 

75.1875 

0520sllb 

19 

75.6570 

0520s11c 

20 

76.4587 

0520bp8j 

21 

75.3466 

0526s2b 

22 

75.5733 

0526s3b 

23 

75.4433 

0603bs5 

24 

76.2651 

0708bs9_ 

25 

75.6427 

0404slb 

26 

75.8893 

0901bs7 

27 

75.7288 

0901bs7b 

28 

75.6417 

0901s7b 
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APPENDIX E. GRAPE AND RVCQ3D INPUT AND PCP CODE 


“GRAPE” INPUT 

&GRID1 

JMAX=340,KMAX=49,NTETYP=3,NAIRF=5,JAIRF=343,NIBDST=7, 

DSI=0.0002,JTEBOT=80,JTETOP=261,XLE=0.0,NOBSHP=7,XTE=1.0, 

ALAMF=0.0,ALAMR=0.0,XLEFT=-.5,XRIGHT=3., 

RCORN=10,NOUT=4,NORDA=5,3,MAXITA=0,300, 

&END 

&GRID2 

AAAI=. 3 ,BBBI=.3 ,DSOBI=0.01 ,ROT ANG=-16.3, 

XTFRAC= 1.5,PITCH= 1.1976,DSRA=.4920,DSLE=.0008,DSTE=.001, 
NLE=48,NTE=36,WAKEP= 1 .,OMEGR=l ,0,OMEGS= 1,0,OMEGP= 1,0,OMEGQ=1.0, 
&END 

&GRID3 AIRFX= 


127.0533, 

127.0025, 

126.9898, 

126.9467, 126.8984, 126.8400, 

126.7714, 

126.6927, 

126.5961, 

126.4768, 126.3752, 126.2990, 

126.2228, 

126.1466, 

125.9942, 

125.7656, 125.3084, 125.0798, 

124.9274, 

124.6988, 

124.5464, 

124.3940, 124.2416, 124.0892, 

124.0130, 

123.9368, 

123.8606, 

123.7844, 123.7082, 123.6320, 

123.5558, 

123.4796, 

123.3906, 

123.0276, 122.6646, 122.3016, 

121.9386, 

121.5757, 

121.2127, 

120.6435, 120.0744, 119.5052, 

118.9361, 

118.3669, 

117.7978, 

117.0168, 116.2359, 115.4549, 

114.6740, 

113.8930, 

113.1121, 

112.1089, 111.1057, 110.1025, 

109.0993, 

108.0961, 

107.0929, 

105.8614, 104.6300, 103.3985, 

102.1671, 

100.9356, 

99.7041, 

98.4585, 97.2128, 95.9671, 

94.7214, 

93.4757, 

92.2300, 

90.9738, 89.7175, 

88.4612, 

87.2049, 

85.9486, 

84.6923, 

83.4351, 82.1778, 

80.9205, 

79.6632, 

78.4059, 

77.1486, 

75.8955, 74.6423, 

73.3891, 

72.1360, 

70.8828, 

69.6296, 

68.3811, 67.1327, 

65.8842, 

64.6357, 

63.3872, 

62.1387, 

60.8997, 59.6607, 

58.4217, 

57.1827, 

55.9437, 

54.7047, 

53.4707, 52.2368, 

51.0028, 

49.7689, 

48.5350, 

47.3010, 

46.0678, 44.8346, 

43.6014, 

42.3682, 

41.1350, 

39.9019, 

38.6769, 37.4520, 

36.2270, 

35.0021, 

33.7771, 

32.5522, 

31.3293, 30.1064, 

28.8835, 

27.6606, 

26.4377, 

25.2148, 

24.0073, 22.7998, 

21.5923, 

20.3848, 

19.1774, 

17.9699, 

16.9731, 15.9764, 

14.9797, 

13.9829, 

12.9862, 

11.9894, 

11.2059, 10.4225, 

9.6390, 

8.8555, 

8.0720, 

7.2885, f 

>.7058, 6.1230, 

5.5402, 

4.9574, 

4.3747, 

3.7919, 3 

1.4079, 3.0239, 

2.6399, 

2.2560, 

1.8720, 

1.4880, 1 

.4199, 1.3437, 

1.2675, 

1.1913, 

1.1151, 

1.0389, 0.9627, 0.8865, 

0.8103, 

0.7341, 

0.6579, 

0.5817, 0.5055, 0.4293, 

0.3531, 

0.2921, 

0.1956, 

0.1270, 0.0762, 0.0381, 

0.0102, 

-0.0051, 

-0.0127, 

-0.0152, -0.0152, -0.0127, 

-0.0051, 

0.0076, 

0.0279, 

0.0559, C 

1.0889, 0.1270, 

0.1702, 

0.2210, 

0.2794, 

0.3454, C 

1.4221, 0.7466, 

1.0712, 

1.3957, 

1.7202, 

2.0448, 2 

1.3693, 2.8953, 

3.4214, 
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3.9474, 

4.4734, 

4.9995, 

5.5255, 6 

.2765, 7.i 

0276, 

7.7786, 

8.5297, 

9.2807, 10.0317, 11.0113, 11.9910, 

12.9706, 

13.9503, 

14.9299, 

15.9096, 

17.1240, 

18.3385, 

19.5530, 

20.7675, 

21.9820, 

23.1965, 

24.4399, 

25.6833, 

26.9267, 

28.1700, 

29.4134, 

30.6568, 

31.9316, 

33.2064, 

34.4812, 

35.7560, 

37.0308, 

38.3056, 

39.6119, 

40.9182, 

42.2245, 

43.5307, 

44.8370, 

46.1433, 

47.4763, 

48.8092, 

50.1422, 

51.4752, 

52.8081, 

54.1411, 

55.4884, 

56.8358, 

58.1831, 

59.5304, 

60.8777, 

62.2250, 

63.5654, 

64.9057, 

66.2461, 

67.5864, 

68.9268, 

70.2671, 

71.5759, 

72.8847, 

74.1935, 

75.5023, 

76.8111, 

78.1199, 

79.3932, 

80.6665, 

81.9398, 

83.2130, 

84.4863, 

85.7596, 

87.0079, 

88.2561, 

89.5044, 

90.7527, 

92.0010, 

93.2493, 

94.4762, 

95.7031, 

96.9301, 

98.1570, 

99.3840, 

100.6109, 

101.8351, 

103.0593, 

104.2835, 

105.5077, 

106.7319, 

107.9562, 

108.9738, 

, 109.9914, 

111.0091, 

112.0267, 

113.0444, 

114.0620, 

114.8767, 

115.6915, 

116.5062, 

117.3210, 

118.1357, 

118.9505, 

119.5655, 

120.1806, 

120.7957, 

121.4108, 

122.0259, 

122.6410, 

123.0433, 

123.4456, 

123.8480, 

124.2503, 

124.6526, 

125.0549, 

125.2271, 

125.3795, 

125.5217, 

125.6487, 

125.7681, 

125.8799, 

125.9815, 

126.0805, 

126.1720, 

126.2583, 

126.3396, 

126.4133, 

126.4844, 

126.5504, 

126.6114, 

126.6673, 

126.7231, 

126.7714, 

126.8197, 

126.8628, 

126.9035, 

126.9390, 

126.9721, 

127.0025, 

127.0279, 

127.0508, 

127.0711, 

127.0864, 

127.1118, 

127.1245, 

127.1118, 

127.0787, 


127.0533, 

AIRFY= 


1.0617, 

0.9855, 

0.9093, 

0.8331, 

0.7569, 

0.6807, 

0.6045, 

0.5283, 

0.4521, 

0.3759, 

0.3200, 

0.2819, 

0.2489, 

0.2210, 

0.1702, 

0.1194, 

0.0838, 

0.0914, 

0.1067, 

0.1448, 

0.1803, 

0.2210, 

0.2667, 

0.3200, 

0.3505, 

0.3810, 

0.4140, 

0.4470, 

0.4826, 

0.5232, 

0.5613, 

0.6045, 

0.6577, 

0.8979, 

1.1299, 

1.3544, 

1.5718, 

1.7825, 

1.9871, 

2.2968, 

2.5944, 

2.8820, 

3.1614, 

3.4345, 

3.7031, 

4.0671, 

4.4257, 

4.7782, 

5.1240, 

5.4623, 

5.7926, 

6.2042, 

6.6019, 

6.9862, 

7.3577, 

7.7168, 

8.0641, 

8.4753, 

8.8712, 

9.2538, 

9.6248, 

9.9861, 

10.3395, 

10.6901, 

11.0331, 

11.3669, 

11.6900, 

12.0012, 

12.2988, 

12.5843, 

12.8558, 

13.1147, 

13.3621, 

13.5993, 

13.8275, 

14.0479, 

14.2594, 

14.4611, 

14.6519, 

14.8307, 

14.9963, 

15.1474, 

15.2840, 

15.4055, 

15.5113, 

15.6010, 

15.6742, 

15.7302, 

15.7704, 

15.7957, 

15.8072, 

15.8060, 

15.7932, 

15.7698, 

15.7350, 

15.6877, 

15.6268, 

15.5514, 

15.4603, 

15.3532, 

15.2293, 

15.0884, 

14.9302, 

14.7545, 

14.5611, 

14.3502, 

14.1228, 

13.8804, 

13.6245, 

13.3564, 

13.0778, 

12.7917, 

12.4966, 

12.1926, 

11.8797, 

11.5580, 

11.2275, 

10.8891, 

10.5426, 

10.1886, 

9.8278, 

9.4608, 

9.0881, 

8.7149, 

8.3365, 

7.9522, 

7.5617, 

7.1646, 

6.7604, 

6.4211, 

6.0769, 

5.7279, 

5.3741, 

5.0158, 

4.6531, 

4.3647, 

4.0722, 

3.7742, 

3.4694, 

3.1562, 

2.8332, 

2.5860, 

2.3330, 

2.0747, 

1.8115, 

1.5438, 

1.2721, 

1.0910, 

0.9086, 

0.7247, 
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0.5397, 

0.3536, 

0.1666, 

0.0457, 

0.0279, 

0.0127, 

-0.0025, 

0.0025, 

0.0178, 

0.1346, 

0.2108, 

0.2870, 

0.5918, 

0.6680, 

0.7391, 

0.9728, 

1.0490, 

1.1252, 

1.4300, 

1.5062, 

1.5824, 

2.5977, 

2.8950, 

3.1850, 

4.7579, 

5.1599, 

5.5503, 

7.4587, 

7.9337, 

8.3922, 

10.4250, 

10.9137, 

11.3861, 

13.4531, 

13.9565, 

14.4440, 

16.2726, 

16.6967, 

17.1081, 

18.6653, 

19.0255, 

19.3718, 

20.6277, 

20.9017, 

21.1577, 

22.0006, 

22.1587, 

22.2936, 

22.5896, 

22.6023, 

22.5911, 

22.3138, 

22.1843, 

22.0292, 

21.1691, 

20.8976, 

20.6055, 

19.2894, 

18.9281, 

18.5564, 

17.0114, 

16.6061, 

16.1909, 

14.4552, 

14.0081, 

13.5599, 

11.7911, 

11.3519, 

10.9113, 

9.3460, 

8.9668, 

8.5851, 

7.2734, 

6.9654, 

6.6594, 

5.6775, 

5.4485, 

5.2154, 

4.4871, 

4.3159, 

4.1398, 

3.7287, 

3.6525, 

3.5763, 

3.2715, 

3.1953, 

3.1191, 

2.8143, 

2.7381, 

2.6619, 

2.3571, 

2.2809, 

2.2047, 

1.8999, 

1.8237, 

1.6713, 

1.0617, 




&END 


0.1321, 

0.0991, 

0.0711, 

0 . 0000 , 

-0.0025, 

-0.0051, 

0.0381, 

0.0660, 

0.0991, 

0.3632, 

0.4394, 

0.5156, 

0.7442, 

0.8204, 

0.8966, 

1.2014, 

1.2776, 

1.3538, 

1.6596, 

1.9804, 

2.2930, 

3.4682, 

3.9135, 

4.3429, 

5.9304, 

6.4574, 

6.9667, 

8.8345, 

9.3880, 

9.9173, 

11.8448, 

12.3975, 

12.9336, 

14.9158, 

15.3830, 

15.8350, 

17.5077, 

17.9059, 

18.2919, 

19.7037, 

20.0282, 

20.3364, 

21.3953, 

21.6180, 

21.8200, 

22.4043, 

22.4910, 

22.5527, 

22.5566, 

22.4993, 

22.4186, 

21.8481, 

21.6457, 

21.4189, 

20.2942, 

19.9741, 

19.6386, 

18.1759, 

17.7956, 

17.4077, 

15.7652, 

15.3362, 

14.8987, 

13.1131, 

12.6704, 

12.2301, 

10.4680, 

10.0967, 

9.7226, 

8.2011, 

7.8922, 

7.5826, 

6.3565, 

6.1300, 

5.9042, 

4.9770, 

4.8173, 

4.6541, 

3.9587, 

3.8811, 

3.8049, 

3.5001, 

3.4239, 

3.3477, 

3.0429, 

2.9667, 

2.8905, 

2.5857, 

2.5095, 

2.4333, 

2.1285, 

2.0523, 

1.9761, 

1.4808, 

1.2903, 

1.1379, 
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“RVCQ3D” INPUT 


'GELDER CONTROLLED DIFFUSION BLADE' 

&NL1 M=340,N=49,MTL=80,MIL= 143 &end 

&NL2 NSTG=4,IVDT=T,IRS=1,EPI=.30,EPJ=.40,CFL=4.0 AVISC2=0., EPS=0.4, 
AVISC4=1. &END 

&NL3 IBCIN=1,IBCEX=1,ITMAX=5000,IRESTI=0,IRESTO=1,IRES=1,ICRNT=10000, 
IXRM=0 &END 

&NL4 AMLE=0.22,ALLE=36.5,BETE=0,PRAT=0.9767,P0IN=1.0,T0IN=L0,G=1.40, 
&END 

&NL5ILT=4,JEDGE=30,RENR=1.0e6,PRNR=.70,TW=0.0,VISPWR= 666666, 
CMUTM=14.0,ITUR=2 &END 

&NL6 OMEGA=0.0,NBLADE=1,NMN=2 &end 

&NL7 TINTENS=0.02 &end 
-.3000 3.000 
1.0000 1.000 
1.0000 .9709 
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“PCP” CODE 


c-c 

C PROGRAM TO READ THE OUTPUT FROM RVCQ3D.F AND GRAPE.F C 
C AND GENERATE A DATA FILE FOR A P/PO VS CHORD PLOT C 
C AROUND THE AIRFOIL. AUG.22, 1991 C 

C-C 

c IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION Q(400,400,4),X(400,400),Y(400,400) 

DIMENSION U(400,400),V(400,400),P(400,400) 

C 

C CALCULATE THE CRITICAL VELOCITY 
cepe=1.00 
rgas=1.0 
pi=1.0 
ti=1.0 

PRAT=0.685 

rhoi=pi/(rgas*ti) 

ceve=cepe-rgas 

g=1.40 

gP=g+l 

cstar=sqrt(2*g*pi/(gp*rhoi)) 

C 

ISTART=80 

IFINIT=261 

C 

READ(7,*)NI,NJ 
print *,ni,nj 

RE AD(7,*)((X(I, J),I=1 ,NI), J= 1 ,NJ), 

& ((Y(I,J),I=1,NI),J=1,NJ) 

READ(7,*)MTL,MIL 

C 

READ(3,*)NI,NJ 

READ(3,*)FSMACH,ALF,RE,TIME 
READ(3,*)(((Q(I, J,K),I= 1 ,NI), J=1 ,NJ),K= 1,4) 

C 

GAMMA=1.4 

C 

DO 1234 1=1 ,NI 
DO 4321 J=1,NJ 

U(I,J)=Q(I,J,2)/Q(I,J,1) 

V(I, J)=Q(I, J,3)/Q(I, J, 1) 

P(I,J)=(GAMMA-1)*(Q(I,J,4)-.5*Q(I,J,1)*(U(I,J)**2+V(I,J)**2)) 

& *RHOI*CSTAR**2 
4321 CONTINUE 
1234 CONTINUE 
C 

C DETERMINE THE XMIN AND XMAX GRID POINT POSITION 
XMIN=0.0 
XMAX=0.0 
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DO 2 I=ISTART,IFINIT 
IF(X(1,1).LE.XMIN)THEN 
IMIN=I 
PRINT *,IMIN 
XMIN=X(1,1) 

END IF 

IF(X(I,1).GE.XMAX)THEN 
IMAX=I 

c PRINT *,IMAX 
XMAX=X(I,1) 

END IF 
2 CONTINUE 

FREE STREAM STATIC PRESSURE USES THE XMIN GRID POINT 
P0=P(IMIN,NJ)+.5*Q(IMIN,NJ,1)*(U(IMIN,NJ)**2+V(IMIN,NJ)**2) 

& *RH0I*CSTAR**2 

i=imin 
j=nj 

t0=l/rgas*((g-l)*q(ij,4)/q(ij,l)-(g-l)**2/(2*g)* 

& (q(ij,2)**2+q(ij,3)**2)/q(ij,l)**2) 
t0=t0*cstar**2 

print *,pO,tO,p(imin,nj),q(imin,nj,l),u(imin,nj),v(imin,nj) 

CALCULATE THE MASS FLOW RATE FROM THE INPUT DATA 
NORMALIZED WITH THE INLET AREA 
dmass=pi/(rgas*ti)*(l/prat)**(-l/g)*sqrt(g*(g-l)*rgas*ti/2* 

& (l-l/(l/prat)**((g-l)/g))) 

print *,'cal. inlet mass flow rate Amass,'* area' 


dmasse=q(imin,nj,l)*RHOI*u(imin,nj) 

print *,'comp. inlet mass flow rate / area =',dmasse 

PRINT OUT DOWN STREAM CONDITIONS 
pl=p(l,l)+.5*q(l,l,l)*(u(l,l)**2+v(l,l)**2)*RHOI*CSTAR**2 
tl = l/rgas*((g-l)*q(l,l,4)/q(l,I,l)-(g"l)**2/(2*g)* 

& (q(l,l,2)**2+q(l,l,3)**2)/q(l,l,l)**2) 

tl=tl*CSTAR**2 
print *,'downstream condition' 
print *,pLtl,p(l,l),q(l,l,l),u(I,l),v(l,l) 

calculate the mass flow rate at exit 
smass=0.0 
ra=0.0 
ba=0.0 

do 21 j=l,nj-l 

smass=smass+(q( 1 j, 1 )+q( 1 j+ L1 ))*(u( 1 j)+u( 1, J+1 ))* 

& RHOI*CST AR*(abs(y( 1 j+ l)-y( 1 j)))*0.25 
ra=ra+(q(lj,l)+q(Ij+l,l))*(abs(y(lj+l)-y(lj)))*0.5* 

& RHOI 

ba=ba+abs(y( 1 j+1 )-y( 1 j)) 
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21 continue 

do 22 j=l,nj-l 

smass=smass+(q(ni j, l)+q(ni j+1, l))*(u(ni j)+u(ni,J+1 ))* 

& RHOI*CSTAR*(abs(y(nij+l)-y(nij)))*0.25 
ra=ra+(q(n 1 j, 1 )+q(ni j+1,1 ))*(abs(y(ni j+l)-y(ni j)))*0.5 
& *RHOI 

ba=ba+abs(y(nij+l)-y(nij)) 

22 continue 

average velocity 
va=smass/ra 
dr=ra^a 

print *,Vel. downsteam,vd =',va 
print *,'ave. density,dr - ,dr 

CALCULATE THE INPUT INLET VELOCITY 
NORMALIZED BY CSTAR 
cO=sqrt(g*rgas*ti) 

vin=sqrt(g*(g-l)*rgas*ti/2*(l-(l/prat)**((l-g)/g)))/c0 
print *,'input inlet velocity / cO - ,vin 

CALCULATE THE CORRECTED DOWNSTREAM VELOCITY 
do 23 j=nj, 1,-1 
vac=u(nij)*vin/va 
alpha=atan(abs(v(nij)/u(nij))) 
vtc=vac*tan(alpha) 
vc=sqrt(vac**2+vtc**2) 
yy=abs(y(ni,nj)-y( 1 ,nj)) 
write(24,*)abs(y(nij)-y(ni,nj))/yy,vc 

23 continue 

do 25 j=2,nj 
vac=u(lj)*vin/va 
alpha=atan(abs(v( 1 j)/u( 1 j))) 
vtc=vac"'tan(alpha) 
vc=sqrt(vac**2+vtc**2) 
write(24,*)abs(y( 1 j)-y(ni,nj))/yy,vc 
25 continue 


CHORD=ABS(X(IM AX, 1 )-X(IMIN, 1)) 

C 

DO 1 I=ISTART,IFINIT 
DIST=ABS(X(1,1)-X(IMIN, 1)) 

XS=DIST/CHORD 

CP=1.0+(-P0+P(I,l))/(.5*Q(IMIN,NJ,l)*(U(IMIN,NJ)**2 
& +V(IMIN,NJ)**2)*RHOI*CSTAR**2) 

wcr=sqrt(u(i,l)' , ‘*2+v(i,l)**2) 

C THE SURFACE POINT P(1,1) ALWAYS EQUAL 0 ? 

PT=P(I,1)/P0 

c WRITE(63,*)XS,WCR 
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c WRITE(64,*)XS,PT 
WRITE(65,*)XS,CP 
1 CONTINUE 
C 

STOP 

END 
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